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Abstract. An exact model for magnetized and rotating outflows, underpressured at their axis, is analysed by 
means of a nonlinear separation of the variables in the two-dimensional governing magnetohydrodynamic (MHD) 
equations for axisymmetric plasmas. The outflow starts subsonically and subAlfvenically from the central grav- 
itating source and its surrounding accretion disk and after crossing the MHD critical points, high values of the 
Alfven Mach number may be reached. Three broad types of solutions are found: (a) collimated jet-type outflows 
from efficient magnetic rotators where the outflow is confined by the magnetic hoop stress; (b) collimated outflows 
from inefficient magnetic rotators where the outflow is cylindrically confined by thermal pressure gradients; and 
(c) radially expanding wind-type outflows analogous to the solar wind. In most of the cases examined cylindrically 
collimated (jet-type) outflows are naturally emerging with thermal and magnetic effects competing in the accel- 
eration and the confinement of the jet. The interplay of all MHD volumetric forces in accelerating and confining 
the jet is displayed along all its length and for several parameters. The solutions may be used for a physical 
understanding of astrophysical outflows, such as those associated with young stellar objects, planetary nebulae, 
extragalactic jets, etc. 
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1. Introduction 

Rotating magnetized objects with hot coronae play a cru- 
cial role in accelerating and collimating astrophysical out- 
flows on the solar, galactic and extragalactic scales. A 
magnetohydrodynamical (MHD) approach is usually con- 
sidered as the first and zeroth-order step for modeling such 
astrophysical outflows. However, the mathematical diffi- 
culties in the treatment of the MHD equations are at such 
a high level that several approximations are still unavoid- 
able to obtain solutions useful for an interpretation of the 
corresponding observed phenomena. 

Several efforts have been put forward recently in 
developing interesting and useful numerical simulations 
(e.g. Sakurai 1985, Ouyed & Pudritz 1997, Bogovalov & 
Tsinganos 1999, Krasnopolsky et al. 1999, Koide et al. 
2000, Keppens & Goedbloed 2000, Usmanov et al. 2000, 
Bogovalov & Tsinganos 2001). Nevertheless, most of them 
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usually apply only to one category of object, such as the 
solar wind, jets from accretion disks, etc, with a limited 
set of explored parameters and boundary conditions be- 
cause the simulations are rather time-consuming. Some 
of these simulations also do not necessarily produce com- 
plete solutions from the base of the outflow up to infinity. 
Others, describe only very transient features because they 
fail to converge towards a stationary or quasi-stationary 
state. And, there are still some doubts on how the bound- 
aries of the numerical box may influence some of those 
simulations, especially in disk winds, a fact related to the 
problem of the correct crossing of the critical surfaces that 
govern MHD flows (Vlahakis et al. 2000). Hence, analyti- 
cal or semi-analytical investigations should be performed 
parallely to numerical simulations, in order to explore a 
more extended set of parameters, despite the necessary 
additional assumptions and simplifications. 

In the framework of steadiness and axisymmetric ge- 
ometry, the MHD system reduces to the well known 
Bernoulli and transfield equations, whose detailed analyti- 
cal treatment is at a rather prohibitive level. General prop- 
erties of these two equations can be studied analytically 
only in the asymptotic regime. For example, in Heyvaerts 
& Norman (1989) it has been shown that collimation into 
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cylinders is a natural configuration for a high speed, su- 
perAlfvenic outflow, in the limit where the field is asymp- 
totically force free but carries a net electric current. A 
complete solution connecting the base of the outflow to 
its asymptotic zone, crossing all singularities, can be ob- 
tained analytically only by neglecting the transfield equa- 
tion and by integrating the ordinary differential Bernoulli 
equation along a streamline. This allows the modeling of 
a 1-D flow around specific regions, e.g., on the equatorial 
plane of a star (Weber & Davis f 967, Belcher & McGregor 
1976), or, along the stellar rotational axis provided that 
the physical variables have been 'a priori' averaged across 
the flux tube (Lovelace et al. 1991). 

An alternative approach largely followed in the last 
years for a 2-D modeling of axisymmetric astrophysical 
outflows is the well known self-similar technique. The ba- 
sis of this treatment is the assumption of a non-linear 
separation of the natural variables of the system of the 
MHD equations. This effectively leads to a scaling law of 
one of the variables as function of one of the coordinates. 
The particular choice of the scaling variable depends on 
the specific astrophysical problem. General properties of 
self-similar MHD solutions are discussed in Vlahakis & 
Tsinganos (1998, henceforth VT98; see also Tsinganos et 
al. 1996). 

Since the well known paper of Blandford & Payne 
(1982), solutions self-similar in the radial direction have 
been investigated to analyse the structure of winds from 
accretion disks (Ostriker 1997, Lery et al. 1999, Casse 
& Ferreira 2000b, Vlahakis et al. 2000, and references 
therein, Aburihan et al. 2001). In most of these models, 
not valid along the polar axis, the driving and collima- 
tion processes derive from a combination of the magnetic 
and centrifugal forces. However, the presence of a hot disk 
corona can also help to drive the outflow very efficiently, a 
factor which may change drastically the initial launching 
conditions (Casse & Ferreira 2000a). 

In a series of studies we have analysed a class of MHD 
solutions that are self-similar in the meridional direction 
(Sauty & Tsinganos 1994, henceforth ST94, Trussoni et 
al. 1997, henceforth TTS97, Sauty et al. 1999, henceforth 
STT99, and references therein). Such a treatment is com- 
plementary with respect to the corresponding radially self- 
similar one because it allows to study the physical prop- 
erties of outflows close to their rotational axis. As in this 
region the contribution to the acceleration by the magne- 
tocentrifugal forces is small, the effect of a thermal driving 
force is very important. Within this model, we either pre- 
scribe the meridional structure of the streamlines, or we 
assume a relationship between the radial and longitudinal 
components of the gas pressure gradient. The main prop- 
erties of the first class of solutions which are asymptoti- 
cally collimated are outlined in TTS97 where the essential 
role of rotation in getting cylindrical collimation is shown. 

If the two components of the pressure gradient are re- 
lated, the meridional structure of the streamlines is de- 
duced from the solution of the MHD equations. Rotating 
magnetized outflows with a spherically symmetric struc- 



ture of the gas pressure may be asymptotically super- 
Alfvenic with radial or collimated fieldlines, depending on 
the efficiency of the magnetic rotator (ST94). This allows 
to deduce a criterion for selecting spherically expanding 
winds from cylindrically collimated jets. In STT99 we ex- 
tended the results of ST94 by performing an asymptotic 
analysis of the meridional self-similar equations assuming 
a non spherically symmetric structure for the pressure. 
It was pointed out that a superAlfvenic outflow may en- 
counter different asymptotic conditions: it can be ther- 
mally or magnetically confined, and thermally or centrifu- 
gally supported. Thus in the frame of this model, the con- 
clusion stated in Heyvaerts & Norman (1989) has been 
extended to cases which are not asymptotically force- 
free, but also include pressure gradients and centrifugal 
terms. Current carrying flows are cylindrically collimated 
(around the polar axis, not necessarily everywhere in the 
flow) for under-pressured flows. While overpressured flows 
attain asymptotically both radial and cylindrical shapes. 

We complete here this study by seeking complete so- 
lutions that connect the basis of the flow with its super- 
Alfvenic regime. In particular we investigate if, and under 
which conditions, the basal region can be matched to the 
asymptotic solutions outlined in STT99. As this analysis 
implies a careful topological study of the MHD self-similar 
equations, with a proper treatment of the critical points, in 
this article we only discuss the structure of underpressured 
outflows (k > 0, see Eq. [?]) that may be asymptotically 
confined by an external pressure or a toroidal magnetic 
field. The study of solutions for pressure supported jets 
(k < 0) is postponed to a forthcoming article. 

In Sec. 2 we outline the main properties of the merid- 
ionally self-similar MHD equations while their asymptotic 
properties analysed in STT9 are summarized in Sec. 3. In 
Sees. 4 and 5 we discuss the main features of the numer- 
ical solutions and their physical validity, while in Sec. 6 
we analyse the dynamical properties of the outflow. The 
general astrophysical implications of our investigation are 
discussed in Sec. 7. 

2. Governing equations for meridional self-similar 
outflows 

We summarize here the main assumptions in a mcridion- 
ally (9—) self-similar treatment of the MHD equations. 
More details may be found in STT94 and STT99, while a 
brief discussion of the various classes of self-similar MHD 
solutions is presented in Appendix A. 

2.1. General properties of axisymmetric steady flows 

The basic equations governing plasma outflows in the 
framework of an ideal MHD treatment are the momentum, 
mass and magnetic flux conservation, the frozen-in law for 
infinite conductivity and the first law of thermodynamics. 
In steady and axisymmetric conditions (Tsinganos 1982, 
Heyvaerts & Norman 1989), the MHD equations allow the 
following definition of the magnetic field B and velocity 
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V through the poloidal magnetic flux function A(r, 9) and 
the poloidal stream function ^(r,9), in spherical coordi- 
nates (r, 0, tp) 
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VA 
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where ^ a = dW/dA L{A) is the total angular momentum 
carried by the flow along a fieldline A = const, fl(A) is the 
angular velocity of this fieldline at the base of the flow, and 



M 



v 2 * 2 
r B 2 A-Kp 



(3) 



is the square of the poloidal Alfven number. From Eqs. ([j]) 
and (||) we see that for trans Alfvenic outflows L(A) and 
Q(A) are not independent integrals: for M = 1, we must 
have L(A)/n(A) = [r 2 sin 2 0] M =i- 

Taking into account Eqs. (Q) and (|2|), the original sys- 
tem of MHD equations reduces to two coupled partial dif- 
ferential equations for the density p and the flux function 
A. 



2.2. Self-similarity: scaling laws of the present model 

There are two key assumptions in a mcridionally self- 
similar treatment of the MHD equations. 

The first is that the Alfven number is solely a function 
of the radial distance, M = M(r). It is then convenient 
to normalize all quantities on the Alfven surface along 
the rotation axis, r = r*. Then denote the dimensionless 
radial distance by R — r/r*, while B», V* and p* are 
the poloidal magnetic field, velocity and density along the 
polar axis at the radius r*, with V 2 = B 2 j>±-Kp*. Also, 
define the dimensionless magnetic flux function a(R,9) = 
2A(r,9)/rlB*. 

The second is the assumption of a dipolar dependence 
of the magnetic flux function A: 



A(r,9)=A(a), 



R 2 



G 2 {R) 



sin 



(4) 



where G 2 (R) is the cross sectional area of a flux tube per- 
pendicularly to the symmetry axis, in units of the corre- 
sponding area at the Alfven distance. Then, for a smooth 
crossing of the Alfvenic surface the regularity condition 
for the toroidal components of the magnetic and velocity 
fields becomes L(A)/Q(A) — r 2 a. 

The next step of the method is to choose the two free 
functions of a (^a and fi) such that the variables (R, 9) 
separate. 

Ratio of magnetic to mass flux. We assume a linear de- 
pendence of ^ 2 (a), which also fixes the density profile, 
Eq. (|): 



^\ = 4717?* (1 + 5a) , 



p(R, a) 



M 2 (R) 



(l+Sa). (5) 



The parameter 6 governs the non spherically symmetric 
distribution of the density. This means to assume a linear 
increase (or decrease) of the density when receding from 
the rotational axis. 

Total angular momentum and corotational speed . For the 
poloidal current within a surface labeled by A(a) = 
const [cx L(a)^ a(o:)] we choose a linear dependence on a 
through the parameter A. Then from the regularity con- 
ditions in Eqs. (|l|) and (0), and the form of '^(a) we 
get 



L(a) — Ar*T4 



il(a) = A 



V* 



1 



r* \/l + Set 



(6) 



Note that A is related to the rotation of the poloidal 
streamlines at R = 1. 

Pressure of the gas . Consistently with the scaling of the 
density, the pressure of the gas is assumed to vary linearly 
with a through the constant k: 



p{R,a) = i P *K 2 n(i?)[i + H 



(7) 



For k > (k < 0) the gas pressure increases (decreases) 
by moving away from the polar axis. We remark that this 
assumption substitutes the usual polytropic relationship. 
Eq. (|) is effectively equivalent with a relationship be- 
tween the heating rate q in the gas, P and A. A special 
case of such a functional relationship among q, p, P and 
A yields the familiar polytropic law P cx p 1 (ST94). 

In this paper we confine our attention to underpres- 
sured jets, i.e., when k > 0. 

Gravity . Consequent to the scaling of the density, the grav- 
itational force per unit volume is 



fa 



GM . 



1 V 2 



M 2 R 2 



(1 + 6a) , 



(8) 



where M is the central gravitating mass and we have 
introduced a new parameter v which is the ratio of the 
escape velocity to the flow velocity on the polar axis at 
R = 1 



GM 

r*V 2 



(9) 



The expansion factor . For homogeneity with the notation 
in ST94 and STT99 we introduce the function F(R), 
which is the logarithmic derivative (with a minus sign) of 
the well known expansion factor used in solar wind theory 
(Kopp & Holzer 1976): 



F(R) = 2 



1 - 



d In G(R) 



dlnR 



(10) 



We remind that the value of F defines the shape of the 
poloidal streamlines (that are parallel to the poloidal mag- 
netic ficldlincs). For F(R) — the streamlines are radial 
(as e.g. in Tsinganos & Trussoni 1991), for F(R) > they 
are deflected to the polar axis (and F — 2 means cylindri- 
cal collimation, see ST94 and TTS97), while for F(R) < 
they flare to the equatorial plane. 
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From the previous definitions we deduce the compo- 
nents of V and B as functions of R, 6, G(R), F(R) and 
M(R). Then, the momentum conservation law provides 
three ordinary differential equations which together with 
Eq. Ji"o| ) can be solved for the four variables M 2 (R), F(R), 
n(i?) and G(R) (see Appendix B). For k > these equa- 
tions have two singularities: on the Alfven surface and 
at the position where the radial component of the flow 
velocity is equal to the radial component of the slow mag- 
netosonic velocity (for a general discussion of the singu- 
larities of the self-similar MHD equations see Tsinganos 
et al. 1996). The details of the regularity conditions that 
must be fulfilled to cross these critical surfaces are widely 
discussed in ST94. They are a crucial element to ensure 
well posed boundary conditions (Sauty et al. 2001). 

2.3. Energetics 

The conserved total energy flux density per unit of mass 
flux density E{a) is equal to the sum of the poloidal ki- 
netic, rotational and gravitational energies, together with 
the enthalpy and heating along a specific streamline. In 
the framework of the present meridionally self-similar 
model, E{a) can be expressed as (ST94): 



(ii) 



The two constants £ and A£ appearing in Eq. (|TT|) rep- 
resent the specific energy along the polar streamline and 
the variation of the specific energy across the streamlines, 
respectively (see STT99). Due to the assumed linear de- 
pendence of the pressure with a, Eq. (0), it turns out that 
the following quantity 



e = AS - k£ 

M 4 



[GRf 



A^ 

G 2 



M 2 - G 
1 - M 2 



M 4 

2N 2 



(8 - K)V 2 



2\ 2 



R 

1-G 2 



1 - M 2 



(12) 



is a constant on all streamlines (ST94). Physically, e is 
related to the variation across the fieldlines of the specific 
energy which is left available to collimate the outflow once 
the thermal content converted into kinetic energy and into 
balancing gravity has been subtracted (STT99). For e > 
collimation is mainly provided by magnetic terms, while 
for e < the outflow is confined mainly by thermal pres- 
sure. Accordingly, in STT99 we defined flows with positive 
or negative e as Efficient or Inefficient Magnetic Rotators, 
respectively (EMR or IMR). 

Note that e is not a mere generalization of the Bernoulli 
constant but rather a variation of the total energy per unit 
mass from one magnetic fieldline to the next. Therefore 
it also contains information on the transfield force bal- 
ance equation and the energies that control the shape of 
the flow tube rather than only information on the various 
driving mechanisms. We also note that an analogous con- 
stant does not exist in the case of meridionally self-similar 



solutions with prescribed streamlines discussed in TTS97, 
where Eq. (^) does not hold. 

3. Asymptotic behaviour of the solutions 

For R 3> 1 the asymptotic parameters of collimated out- 
flows (Foe — 2, Goo and M m bounded) depend on the 
value of e and the force balance across the poloidal stream- 
lines, /vp + /b + fc = 0, with / V p, /b and f c the pres- 
sure gradient, magnetic stress and centrifugal force, re- 
spectively. Then, we may calculate and Goo as func- 
tions of the parameters e/2A 2 , k/2X 2 and Hoo- 

The asymptotic properties of these self-similar winds 
have been discussed in detail in STT99, and here we sum- 
marize their main features for the case of underpressured 
outflows (k > 0). 

— Two main confining regimes exist. In one the outflow 
is collimated by the pinching of the toroidal magnetic 
field (e > 0, magnetic regime: EMR) and in the other 
by the thermal pressure (e < 0, thermal regime: IMR). 

— For k — ► we have collimation only for e > 0: there 
is no pressure gradient across the streamlines and the 
flow can be confined only by the magnetic stress (/#, + 
fc = 0) 

— The flow is supported by the centrifugal force and col- 
limated cither by the thermal pressure (e/2A 2 -C 0, 
fc + /vp = 0), or by the magnetic pinch (e/2A 2 ~ 0, 
fc + fs= 0). 

— The collimated streamlines always show oscillations. 
This behaviour is consistent with the results found in 
more general, non self-similar treatments (VT98). 

In conclusion and within the present model, from the 
asymptotic analysis it turns out that underpressured and 
meridionally self-similar outflows should be in principle al- 
ways collimated. However, for very small values of k and 
asymptotically vanishing pressure — 0, an asymptoti- 
cally radial configuration of the fieldlines is not excluded. 

4. Numerical results 

4.1. Numerical technique and parameters 

Using routines of the NAG scientific package suitable for 
the treatment of stiff systems and the Runge-Kutta algo- 
rithm, Eqs. ([To|), (B.7) - (B.13) are integrated upstream 
and downstream of the vicinity of the Alfven point (Ri n — 
1 ± dR) with M m = 1 ±pdR and G in = 1 ± (2- F in )dR 
( F in ~ i**). The slope p of M at R = 1 is given in Eq. 
( B.14[ ) . We first integrate upstream tuning the value of F in 
until we select the critical solution that smoothly crosses 
the singularity and reaches the base of the wind R , where 
M — ► 0. With this value of F in we then integrate down- 
stream to the asymptotic region (R^ usually between 10 4 
and 10 6 ) and we get the complete solution. The value of 
the pressure n ira (ss IT*) is chosen such that IT(i?) is pos- 
itive everywhere (see however Sec. 5 for a possible release 
of this constraint). 
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We have seen in Sec. |^ that the asymptotic properties 
of collimated outflows are ruled by the three parameters 
e/ 2A 2 , re/2A 2 and ILx;. As it is evident from Eqs. (B.7) 
- ( [B.9D and (|T2|), the numerical solutions require five pa- 
rameters: e, re, A, 6 and v. However only four of them arc 
independent due to the constraint imposed by the integral 
e, Eq. (113) , which has the following expression at R = 1: 



( K _ S)u 2 + A 2 (r 2 + 1) - (1 - re) + F 2 /4 . 



(13) 



where r[= (2 - F*)/p] is given by Eq. flB.14p . 

We remind that the parameters A, 6 and re are related 
to the strength of the rotational velocity, and to the struc- 
ture of the density and of the pressure in the 9 direction, 
respectively. For 6 and re positive the density and pressure 
increase when receding from the polar axis. From the ex- 
pression giving the base of the outflow R a as deduced from 
the definition of e with M = and G < 1, Eq. Q: 



R = 



{s- K y/\ 2 

2 - G 2 - e/A 2 



(14) 



it is evident that the quantity (5 — k)v 2 rules the flow dy- 
namics in the subAlfvenic region. Large values of it lead 
to steep initial acceleration, with the base of the flow close 
to the Alfven surface. If 6 = re, Ro = and there is no ac- 
celeration close to the base for a given speed at the Alfven 
transition. We also see that we must have re < 5 in order 
that Rq > while S > (S < 0) increases (decreases) the 
initial acceleration. For re > 5 the initial thermal driving 
of the outflow is not possible at the base (an analogous 
case was found in TTS97, where the quantity 8v 2 had to 
be larger than a minimum threshold value to have mass 
ejection). 

In the following two subsections we have fixed A = 3 
and 5 = 0.01 [with the value of v deduced from Eq. |l3| 
and analysed the trend of the solutions for different values 
of e/2A 2 , re/2A 2 and II*. A summary of all these properties 
is given in subsection (4.4)forA = ^= l. 



4.2. Behaviour of the solutions with re 

In Figs. |l| we have plotted the Alfven number M(R), the 
radial velocity V r (R) along the polar axis and the jet ra- 
dius G(R) for e/2A 2 = 0.25 and -0.1, and four different 
values of re = 0, 0.001, 0.003, 0.009. The chosen value of 
LT* is the minimum allowed (II* im i n ), i.e. such that II is 
either vanishing asymptotically, or it is zero in the deepest 
minimum of the corresponding oscillations of the stream- 
lines. With these assumptions we deduce the values of 
v w 20 4- 100. The main parameters of the solutions are 
listed in Table 0. 

In the subAlfvenic region (R<1), independently of the 
value of e, the increase of re rises the value of F* (i.e. the 
rate of expansion of the streamlines at R = 1 is lower) and 
of the pressure: the outflow is 'hotter' and 'narrower' in the 
transAlfvenic region. For re = — > 0.009 n* m i n « 200 -j- 
2000 (see Table These effects however are evident only 
for e < 0, while for positive e the solutions are practically 



unaffected by the value of re (unless re — * 5, in which case 
oscillations of small amplitude occur and R decreases). 
The location of the base of the flow R is controlled by 
(<5— n)v 2 . However in our plots the base of the flow does not 
change drastically except when re — * 8 (see Table [l]). These 
trends are basically unaffected by the value of e/2A 2 . 

In the super Alfvenic region (R>1), when e > the 
results are basically insensitive to an increase of the values 
of re. The flow is always collimated as in ST94, while the 
radius of the jet is smoothly increasing. On the other hand, 
when e < 0, the effect of re ^ is much more drastic. 
The outflow rapidly expands with R, reaches a maximum 
width and then it sharply recollimatcs to a rather narrow 
jet with large amplitude oscillations. For very small values 
of re (^cC S) this bump increases its amplitude while its 
position moves downwind. In the limiting case of re — > 0, 
the bump is shifted to R = oo and the flow has radial 
streamlines. For re — > S the bump and the oscillations have 
shifted upstream and smoothed, and the solution becomes 
similar to the one with e > 0. Note also that at the position 
of the bump the Alfven number has a maximum (and the 
density is minimum), while the velocity is monotonically 
increasing and attains asymptotic values larger than in 
solutions with e > 0. 



4.3. Behaviour of the solutions with II* 

The effect on the solutions of increasing the pressure at 
the Alfven distance II* to values II* > II* !m i n is shown 
in Fig. |, where G(R), U(R) and T(R) are plotted for a 
positive and a negative value of e. First, we note that only 
for II* 3> n* m i n the solutions with positive and negative 
e appear to have similar profiles. And second, it is also 
evident that an increase of II* always leads to a reduction 
of the jet radius, i.e. it has basically the same effect as 
the increase of re (Fig. 1). This is expected because the 
pressure gradie nt a cross the streamlines has terms which 
are cx rell, Eq. flC.7| ). Nevertheless, the asymptotic Alfven 
number and density are scarcely affected by the value of 
II*. And, the flow velocity (oc M 2 /G 2 ) increases with II* , 
a fact also related to an increase of the thermal driving 
efficiency. 

The pressure at the Alfven surface H» cannot be how- 
ever arbitrarily large. Namely, the asymptotic pressure 
LToo increases with LT*, and it may even exceed the corre- 
sponding pressure values at R = 1 , but there is a sharp de- 
crease of II at the position of the bump (with a consequent 
lowering of the temperature). This minimum value of the 
pressure at the bump decreases by rising II*, and above 
a maximum value (II* ]max ) the pressure becomes negative 
at the bump. Hence, acceptable solutions exist only for 
n* jm in < II* < n*. max . In other words, the range of ac- 
ceptable values of II* is quite narrow, and further shrinks 
by decreasing e/2A 2 and increasing re or 5 (see Table 1). It 
turns out that below a threshold value of e/2A 2 the solu- 
tions are unphysical because we obtain II* jm i n > n* max , 
unless re — > 5. For example for e/2A 2 = —5 we have accept- 
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Fig. 1. Plots of M (upper panel), V r (in units of V*, middle panel) and G (lower panel) vs R along the polar axis for 
e/2A 2 = -0.1 (left panel) and 0.25 (right panel). Labels refer to different values of k = (a), 0.001 (b), 0.003 (c) and 
0.009 (d). In all these solutions we have selected the solution corresponding to II* = n*, m i n . 



3 and 8 — besides the solution with n = and radially expanding 
streamlines. 

Concerning the asymptotic velocities of the collimated 
outflow, with the assumed parameters, as k = 0.001 — > 
0.009 and II* = n* ;m .in — * n* iinax we find for e/2A 2 = 0.25, 
M = 3.5 + 8 and V r /V* = 5 4- 5 x 10 2 , while for e/2A 2 = 
-0.1 we get M = 3 15 and V r /V* = 20 4 x 10 2 . 



Table 1. Parameters of the solutions for A = 
0.01. R and R cr are given for II* = n* jm i n . 



e/2A 2 = -0.1 



ft 


n*,min 




Ro 


Rcr 


V 





289 




0.72 


0.72 


34.1 


0.001 


321 


505 


0.73 


0.73 


35.9 


0.003 


413 


653 


0.73 


0.74 


40.3 


0.009 


1951 


2089 


0.79 


0.97 


102 


e/2A 2 = 0.25 


K 


n*,min 




Ro 


Rcr 


V 





214 




0.60 


0.62 


24.1 


0.001 


236 


2007 


0.61 


0.62 


25.2 


0.003 


297 


1135 


0.61 


0.63 


28.2 


0.009 


1286 


1423 


0.45 


0.88 


65.1 



able solutions only for k = 0.009 with 4750 < IT, < 4835, 



4.4. Summary of the parametric analysis 

In order to summarize the trends we just discussed, we 
illustrate in Fig. ^ how the behaviour of the solutions with 
the pressure is displayed in the parameter space (e, k) for 
an arbitrary and illustrative set of values A = v = 1. The 
exact location of the various domains in this space [c,k] 
depends strongly on the set of the remaining parameters. 
The following discussion in each case refers to the solution 
corresponding to the minimum pressure II* , m in allowed, 
except where it is otherwise indicated. 
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Fig. 2. Plots of G (upper panel), II (middle panel) and T (lower panel) vs R along the polar axis for k = 0.003 and 
e/2A 2 = —0.1 (left column) and e/2A 2 = 0.25 (right column). Labels refer to different values of II*. For e/2A 2 = —0.1, 
II* « 413 (a, n*, min ), 420 (b), 653 (c, II*, max ). For e/2A 2 = 0.25, IF « 297 (a, n*, min ), 400 (b), 1135 (c, II„, max ). 



In the domain of Efficient Magnetic Rotators (EMR, 
e > 0), jets are cylindrically collimated by magnetic forces, 
as expected (domain with stars in Fig. |J), unless k is get- 
ting too large in which case we enter the regime of pres- 
sure confinement (crosses in Fig. ^). This pressure con- 
fined region extends to the domain of Inefficient Magnetic 
Rotators (IMR) on the left part of Fig. || However, when e 
is getting too negative, it occurs that II* , mm > II* , max for 
small enough values of k (open squares in Fig. ||). This re- 
gion corresponds in principle to physically non acceptable 
collimated solutions; we shall see however in the following 
subsection that a subdomain exists where radial solutions 
may exist. 

We remark that these main features of the numerical 
solutions we just discussed are not qualitatively affected 
by different values of the parameters. The main effect of 
lower values of A is to shift downstream the 'X' type crit- 
ical point closer to the Alfven point (see ST94), and the 



properties of the solutions appear very similar for equal 
values of the quantity (S — k)v 2 . 

In the previous subsections, we have seen that slight 
increase of the magnitude of k and/or II* sharply reduces 
the jet radius (mainly for e < 0). The jet radius may 
even asymptotically become smaller than at R — 1 (i.e., 
Goo < !)• This configuration implies a reversal of the 
electric current somewhere in the superAlfvenic region, 
and this may still be consistent with a reasonable topo- 
logical structure of the outflow (see e.g. VT98, Vlahakis 
1998). However jets with such a decreased transverse size 
far from the central object are likely to be unrealistic. 
Therefore we could assume in general that the parameters 
should be constrained such that Goo > 1. We must note 
also that when Goo < 1 the asymptotic pressure is rather 
high (see Fig. ||) while it is reasonable to assume that in 
astrophysical jets Hoc <C n*. From all these arguments it 
appears that collimated solutions suitable to model colli- 
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(IMR) 



-0.02 0.02 
e / 2A 2 



0.05 0.1 0.2 



(EMR) 



-0.02 0.02 
s / 2A 2 



0.05 0.1 0.2 



Fig. 4. Summary of the various domains of the pres- 
sure behaviour in the parameter space [e/2A 2 , k/2A 2 ] for 
v = X = 1. We show the domains of magnetic collimation 
(stars), pressure confinement (crosses), the forbidden re- 
gion where all collimated solutions obtain at the bumps a 
negative pressure (open squares), and finally the domain 
of the radial solutions (filled squares). 



Fig. 5. In the same parameter space as Fig. 4 we show the 
limiting solutions between the collimated solutions and 
the flaring ones which are either collimated without os- 
cillations but negative IToo (circles) or which are radially 
expanding at least up to R = 10 6 with II^ = (filled 
squares) 



mated astrophysical outflows should be in general selected 

with >n* im i n . 
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n > n 

* *,lim 

n = n 

* *,lim 



n < n 

* *,lim 



<3 



b) 



n < n 



*,lim 



n = n 

* *,lim 



n > n 

* *,lim 



R 



R 



Fig. 6. U(R) (panel a) and G(R) (panel b) are plotted for e/2A 2 = -0.5, k = 0.001, v = A = 1. By tuning II* to a 
special value II* = II*ji m we find, between a sharply recollimated (II* > Il^nm) and a flaring solution (II* < n*ji m ), 
a limiting solution (dashed line) which is radial up to R ~ 10 6 (see the plot of the streamlines in Fig. @a). 




Fig. 7. The same as in Fig. 6 for e = 0, k — 0.01, v = X = 1: the limiting solution now is collimated without bumps. We 
plot also the solution with positive pressure everywhere, labeled with II*, m i n . The streamlines of the limiting solution 
and of the one with II* > IFjim are shown in Figs. |p and He, respectively. 



5. Limiting solutions 

We outline here the properties of two kinds of solutions 
(that we call 'limiting' solutions) that are outside the 
boundary conditions on the pressure we discussed in the 
previous section. 

5.1. Radially expanding limiting solutions 

We have seen that for sufficiently negative values of e and 
small values of k, the sharp jet recollimation leads to a neg- 
ative pressure, with II* min becoming larger than II* max . 



These solutions have been rejected as unphysical (open 
squares in Fig. ^). However it turns out that by further 
decreasing II* below n*, m in we find a subdomain of this 
region (filled squares in Fig. [^), where we could follow nu- 
merically solutions that remain radial up to R = 10 6 by 
tuning the initial pressure to a specific value II* = IFjim. 
An example of such a solution is given in Figs. ^| and ||a. 
These solutions, with an almost vanishing pressure and 
radial asymptotics, are limiting solutions between sharply 
recollimated ones (for II* > H*fi m , see Fig. |^) and flaring 
solutions (for II* < II* i im ). Numerically we are unable 
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to say if such solutions are strictly radial asymptotically, 
or if they become cylindrical with a negative pressure far 
away from the base, as we discuss hereafter. However they 
can model radial winds up to the region of the outer as- 
teropause where the wind is shocked with the interstellar 
medium because we do not expect the shock to be that far 
away (i.e. farther than R = 10 6 ). Thus they allow some 
extension of the domain of radial winds into that of under- 
pressured flows and they are well adapted to model winds 
like the solar wind (Lima et al. 2001a). 

5.2. Non oscillating limiting solutions 

In the remaining region of Fig. 5 (circles) the limiting 
solution between flaring and recollimating ones is a solu- 
tion that collimates into cylinders without oscillations or 
bumps, and with negative n^,. An example is presented 
in Figs. [7] and ||b. This result is not surprising and can be 
related to that obtained by TTS97, where two families of 
solutions were found for fixed streamlines. The first one 
was composed of thermally confined and underpressured 
jets showing the presence of a bump in the pressure. The 
second family was composed of initially underpressured 
jets which become overpressured downstream and mag- 
netically confined without bumps. We have the analogous 
result here, at least in the domain of negative e. 

But can this solution with negative II be acceptable? 
In fact, the pressure can be defined modulus an arbitrary 
constant, 



P(R, a) = \p^IV{R)[l + ku] + P (> 



(15) 



Thus we can drop the constraint II* > IT jm i n by choosing 
values of this constant such that the pressure is always 
positive. However the increase of P implies also higher 
temperatures of the plasma, that should be eventually de- 
tectable with observations. Furthermore, high gas pressure 
could require unreasonable large heating processes in the 
plasma. So adding P Q to solutions that present a bump 
in the pressure (as in Fig. ||c) is likely to be unrealistic. 
Conversely, including this constant to the limiting solu- 
tions with no bump (e.g. Fig. ||b) is feasible in principle 
and should not require a much larger amount of heating. 
We could speculate that such solutions could be interest- 
ing to model jets which are initially underpressured (in 
the region where n > 0) and overpressured in the outer 
zone (where n < 0). 

6. Dynamics of the outflow 

6.1. Force balance along and across the flow 

A complete physical understanding of the behaviour of 
the solutions discussed in the previous section may be ob- 
tained by studying systematically the various forces act- 
ing on the outflow. In STT99 a similar analysis was at- 
tempted in asymptotically cylindrically collimated out- 
flows. In such a case the jet had attained its terminal 



speed and force balance was studied only in the normal 
to the fieldlines direction. In the present case we consider 
all forces across and along the flow, in the whole region, 
from the base to the asymptotic zone. For this purpose, 
the original MHD set of equations can be projected on 
the poloidal plane and there they may be split in two 
directions, one being tangent (s) and the other (n) per- 
pendicular to a particular poloidal streamline a = const. 
In this local system of orthogonal coordinates, the most 
general form of the MHD equations reduces to: 







-pv P 



dVp 
ds 



dP 
~to 



dV 

'to 



Attzu 



dm 
~to 



II n : 



pv p 2 


dP 


dV 
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dn 


on 


+ ( 




d ( 


B2 A 


d 




dn { 


8tt I 


dn 




B 



where w — rsin(9, V is the gravitational potential (= 
— QM.fr) and R c is the radius of curvature of the poloidal 
streamline. The vectors s and h are directed towards the 
outer regions and the polar axis, respectively. Then, a force 
/ has two components, f s and /". A positive component 
f s accelerates the flow along the streamline, while a posi- 
tive value for /" collimates the plasma. Conversely, nega- 
tive f s and /" imply that the wind is decelerated and that 
the streamlines 'flare' away from the axis, respectively. 

The full expressions of the forces are given and dis- 
cussed in detail in Appendix C. In the following we out- 
line the main dynamical properties along and across the 
streamlines, referring to the cases of collimated flows we 
have discussed in Sec. ^ and plotted in Fig. ||. 

6.2. Force balance along s 

It is easy to see that, the centrifugal force and the neg- 
ative gas pressure gradient for a monotonically decreas- 
ing H(R), always accelerate the flow (see left panels of 
Figs. H - |l^). On the other hand, gravity and the hoop 
stress associated with the toroidal magnetic field always 
decelerate the plasma. However the initially negative gra- 
dient of B^/8tt dominates close to the base and acceler- 
ates the flow (see left panels of Figs. ||- |l^). Finally, as is 
well known, the poloidal component of the magnetic field 
has no effect on the plasma dynamics along s. From the 
left panels of Figs. || - [l2], we see that the outflow is al- 
ways thermally driven, with the magnetic and centrifugal 
forces playing a rather marginal role in the acceleration. 
For the assumed parameters the overall acceleration oc- 
curs for R<10 3 , independently of the values of e and k. 
Oscillations in the flow variables may exist in some cases. 
For example, when e < and for intermediate values of 
k (Fig. O), asymptotically the pressure gradient reacts 
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Fig. 8. Sketch the poloidal streamlines of the radial limiting solution of Fig. || (panel a). In panels (b) and (c) we show 
the poloidal streamlines of the limiting solution (collimating without oscillations) and of the one sharply recollimated 
with n» > n^iim, reported in Fig. ^. 



to the oscillations of the streamlines, as expected in ther- 
mally confined jets (STT99). For larger values of k the 
flow appears already collimated in the acceleration region 
(see below), so the oscillations appear much closer to the 
Alfvenic surface (Figs. |l^ and [l2[ left panels). 



6.3. Force balance along h 

The analysis of the forces acting perpendicularly to a 
poloidal streamline is more complicated. First, the cen- 
trifugal force is always negative tending to open the 
streamlines. Second, both the tension and the gradient of 
the toroidal magnetic field are directed towards the axis, 
assisting collimation. Also, the normal component of the 
gravitational field tends to collimate the flow when the 
streamlines do not flare faster than the spherical expan- 
sion. Two terms contribute to the gradient of the pressure, 
with opposite effects: the term related to the radial compo- 
nent (oc dn/di?) decollimates the fieldlines while the one 
related to the meridional component (oc kH) is positive 
and confines the flow. Concerning finally the poloidal mag- 
netic field, its gradient is negative (decollimating) while its 
tension is positive (collimating). 

The interplay of all these forces is shown in the right 
panels of Figs. || - [l^. For positive e and intermediate 
k (Fig. |^) the dynamics of the streamlines in the sub- 
Alfvenic and transAlfvenic region (R < 10) is governed by 
gravity and the radial component of the pressure gradi- 
ent. Farther away the magnetic and centrifugal forces pre- 
vail: the centrifugal force is balanced by the hoop stress 
and the gradient of the toroidal magnetic field (these two 
components provide the same contribution to this force). 
Asymptotically then the flow is magnetically confined (see 
also STT99). The effect of the other forces is quite negli- 



gible: they just show an oscillating behaviour that slightly 
affects the whole structure. This same picture basically 
holds if we increase the value of k, with an increased ef- 
fect of the thermal pressure (Fig. |l^, right panel). Namely, 
the wind is always asymptotically magnetically collimated 
but in an intermediate region (R w 15) the outflow is 
thermally confined. Oscillations are always present in the 
subAlfvenic and the intermediate regions with larger am- 
plitude: this is likely to be related to the larger effect of the 
pressure on a more collimated wind in this intermediate 
region. 

A similar scenario is also found for negative e and large 
values of k (see Fig. (L2L right panel). The wind is also 
asymptotically magnetically collimated but with a larger 
intermediate region (R sa 20 4 100) where the outflow is 
thermally confined. 

For smaller n (Fig. [ll], right panel), differently from 
the case with positive e, the toroidal magnetic force bal- 
ances the curvature force in an intermediate region (R ss 
104 100) but without any confinement of the flow that ex- 
pands almost radially in this zone. The collimation starts 
downwind where the gas pressure prevails on the magnetic 
force (R ss 100 4 400), leading to a cylindrical configu- 
ration with oscillations (see Fig. |^a). In this asymptotic 
region the centrifugal force is balanced by the gradient of 
the pressure, with a small contribution from the toroidal 
magnetic field, since we are there in the thermally confined 
regime. 

7. Summary and astrophysical implications 

In this paper we confined our attention to the study of 
outflows with a density increasing from the axis towards 
the surrounding streamlines [cf. Eq. (||) with S > 0] faster 
than the pressure does [cf. Eq. (Q) with k > 0], i.e. 5 > k. 
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Fig. 9. Plot of the different forces acting along (I, left panel) and perpendicular to the poloidal streamlines (n, right 
panel) for k = 0.003, e/2A 2 — 0.25 and II* = n m i n (Fig. ||b). We plot the negative of the pressure gradient (dashes), 
the gravitational force (dash-three dots), the poloidal inertial and curvature forces (dots), the poloidal magnetic force 
(dash-dot), the centrifugal force (thin solid line) and the total toroidal magnetic force (thick solid line). 









a) 


J v 


'VP 










-rftar .-1 iH.: 






vji i mi — m <m^u^ 























10* 




Fig. 10. The same as Fig. | for n = 0.009, e/2A 2 = 0.25 and II* 



n, 



(Fig. |d). 



In such underpressured outflows, the temperature is 
peaked at the axis relatively to the surrounding regions. 
Such centrally hotter astrophysical atmospheres are likely 
to exist in coronae wherein, for example, strong Alfven 
waves propagating along the direction of the polar mag- 
netic field deposit enough energy which raises thus the 
plasma temperature. The inevitable consequence of this 
high temperature is an induced wind-type outflow, which 
is dominantly "thermally" driven around the rotation axis 
where magnetocentrifugal forces are negligible (see Figs. ^| 
- O ). By thermally we mean not only thermal conduction 



but also waves, etc, which can be effectively included in 
the pressure gradient term. 

After the first acceleration stage where the gas expands 
almost radially up to the transAlfvenic region, if the out- 
flow carries an electrical current and the magnetic rotator 
is sufficiently efficient (EMR), it is collimated magneti- 
cally via the combination of the hoop stress and the pres- 
sure gradient of the azimuthal magnetic field in the super- 
Alfvenic regime. Several jets from Young Stellar Objects 
have probably such smoothly increasing radius with a 
large terminal transverse extension. They would naturally 
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Fig. 11. Plot of the different forces acting along (s, left panel) and perpendicular to the poloidal streamlines (n, right 
panel) for k = 0.003, e/2A 2 = —0.1 and II* = II m i n (Fig. |3^,). As in Fig. ^| we plot the negative of the pressure 
gradient (dashes), the gravitational force (dash-three dots), the poloidal inertial and curvature forces (dots), the 
poloidal magnetic force (dash-dot), the centrifugal force (thin solid line) and the total toroidal magnetic force (thick 
solid line). 
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Fig. 12. The same as Fig. [ll] for k = 0.009, e/2A 2 = -0.1 and II* = n min (Fig. |c). 



correspond to the magnetic collimation produced by an 
EMR. Outflows from an EMR may also describe jets from 
Radio Loud AGNs like FR Is and FR lis (see Sauty et al. 
2001 for a possible application of the present solutions to 
AGNs). For definite conclusions in this last case, however, 
a relativistic extension of our model is needed. 

On the other hand, the produced outflow may also be 
collimated via thermal pressure gradients, if it is under- 
pressured on the axis. Such pressure confined jets tend 
to show a rather strong recollimation. This could be re- 



lated with some jets from slowly rotating T Tauri stars 
and Planetary Nebulae which show a rather strong recol- 
limation, or some choked winds from Seyfert Galaxies. For 
example, the peculiar jet of RY Tauri originates in a very 
slowly rotating CTTS and seems to recollimate at 38 stel- 
lar radii from the base (see Gomez de Castro et al. 2001). 
The sharp gradient of the gas pressure in the recollimat- 
ing region could lead naturally to the formation of a shock 
which may be consistent with the scenario of recollimation 
by internal shocks proposed in the literature. 
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In very Inefficient Magnetic Rotators (IMR) like the 
sun, radial solutions as the one plotted in Fig. |^ are likely 
to apply to winds produced by such sources (Lima et al. 
2001b). Thus, jets from Young Stars may evolve from a 
narrowly collimated magnetic outflow to a radially ex- 
panding wind. Such radial solutions could also correspond 
to winds from Seyfert Galaxies. 

If we were to apply recollimating solutions to jets from 
Planetary Nebulae, we would be inclined to think that in 
these objects recollimation is an effect of the pressure gra- 
dient rather than a pure magnetic pinching. Nevertheless, 
the jet could be ultimately magnetically confined once it 
has refocalized. 

Interesting solutions for underpressured jets becoming 
overpressured after exiting from the central embedding 
medium can be found. They correspond to the limiting 
solution plotted in Figs. and |^b. Such outflows are mag- 
netically confined and in contrast with the previous ones 
do not exhibit any recollimation or oscillations. They may 
describe Young Stellar jets exiting from the central cloud. 

More precise modelling of various specific astrophysi- 
cal objects is underway and will be presented elsewhere; 
here we just summarize some of the results obtained so 
far. Preliminary applications to the solar wind are given 
in Lima et al. (2001b). By analysing Ulysses data we ar- 
rived at the following ranges of the model parameters for 
the solar wind : 5 pa 2 - 4, A w 0.1, « pa 0.3 - 0.8 and 
v pa 0.7. With these parameters a reasonable fitting of all 
observed quantities is obtained, except perhaps the mag- 
netic field at large distances, which turns out to be a bit 
too weak compared to the measured one. The above values 
of the parameters imply e/2A 2 ~ —50, thus confirming the 
strong inefficiency of the solar magnetic rotator. In such a 
case the solutions are expected to be qualitatively similar 
to those shown in Fig. 3a, i.e. radially expanding stream- 
lines connected through the bump to the far dowstream 
collimated region. In the present case, with k <C S and very 
low asymptotic pressure, the position of the bump is very 
far away from the Alfvenic surface, then the streamlines 
are expected to be basically radial in the region explored 
by the spacecraft. 

Parallel to this, preliminary results for jets from T 
Tauri stars (Meliani 2001) provide the following range 
of the values of the parameters: S pa 0.1 — 0.2, A pa 0.8, 
K pa 0.01 - 0.06, v w 1 - 6 and Pa 10 Pa 20 - 30 r . 
It is worth noting that observations of these objects im- 
ply the existence of a UV emitting region near the central 
star, which may possibly be associated with a shock (e.g. 
in RY Tau this UV zone is at pa 38 r*). With the above 
parameters we deduce e/2A 2 pa (negative or positive), 
which means that T Tauri stars correspond to more effi- 
cient magnetic rotators. For the above ranges of the pa- 
rameters, the solutions are expected to be similar to those 
reported in Fig. 7 and 8. In particular we could have a rec- 
ollimated solution with the UV shock at the position of 
the bump, or a limiting solution with monotonically col- 
limated streamlines without oscillations (see Fig. 8b), in 
which case the shock would result from the growth of some 



internal instability. Both types of solution reproduce in a 
satisfactory way the velocity and density profiles, the jet 
morphology and the rotational rate. Yet, the second type 
of solutions leads to lower temperature profiles which sug- 
gests that a thermally driven wind could account for the 
observed data in such objects, even without extra non- 
thermal heating. 

To distinguish between recollimated and limiting so- 
lutions, we may ask about the stability of our solutions. 
Some hints about this crucial point can be given for the 
asymptotic region. By a local analysis using the thin flux 
tube approximation, Hanasz et al. (2000) have shown that 
the internal part of jets with rotational laws comparable to 
the one used in the present model (Eq.^) are more stable 
to magnetorotational instabilities than the surrounding jet 
coming from the disk. In some cases rotation completely 
suppress the instability. Of course, we do not claim that 
this analysis is complete but it may be an indication that 
azimuthal magnetic fields do not necessarily destroy the 
global structure of MHD outflows modelled by our solu- 
tions. 

However, while the properties of the local stability of 
cylindrical jets are quite well known in the linear regime 
and for simple equilibrium structures (see e.g. Birkinshaw 
1991, Ferrari 1998, and references therein, Kim & Ostriker 
2000), unfortunately a global, non-linear analysis of the 
stability of general axisymmetric MHD steady flows is not 
yet available. This still remains a challenge for future stud- 
ies. Furthermore, the present solutions may be used for the 
testing and interpretation of numerical codes employed for 
such a stability analysis. 

To conclude with, let us stress that in the present study 
cylindrical jets and radial winds can carry a net current, 
at least in the central part close to the axis, because they 
are not necessarily force-free asymptotically. This is not 
in contradiction with the main conclusion in Heyvaerts 
& Norman (1989), because some of the approximations 
done there do not necessarily hold close to the axis. Our 
description is also only valid at the central part of the out- 
flow, and the current should close either in the surrounding 
wind, the external medium, or in a current-carrying sheet 
around the jet. 
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Appendix A: Generalities on self-similar solutions 

For the sake of generality the free functions A(a), ^ a (ch) 
and f2(a), can be defined through the following new func- 
tions (for details see VT98): 

gi(a)(x J Fi(A')da, g 2 (a) cx J flH 2 A da, 

93(a) cx T 2 (A.l) 

where A' = dAjda. We remark that g\ is related to the 
magnetic field structure through the magnetic flux func- 
tion, while g 2 and 173 define the current and density distri- 
bution, respectively. Suitable choices of T\ and T 2 allow 
to select various classes of self-similar solutions. 



Radially self-similar models . If we assume M = M(9) (i.e. 
the surfaces with the same Alfven number are conical), 
G= G{9) and in Eqs. (A.l): 



A' 2 

^(A') = — , T 2 {^ A ) = 
a 



A -da . 



,3/2 



(A.2) 



we obtain a class of solutions that are self-similar in 
the radial direction, used to model outflows from disks. 
Assuming for instance gi, g 2 , gs oc a x /(x — 2) (see Table 3 
in VT98) with x = 3/4 we have the well known solutions 
of Blandford & Payne (1982). 

Meridionally self-similar models . By choosing G = G(r), 
M = M[r) (i.e. the surfaces with the same Alfven number 
are spherical) and: 



^(A'j^A* 2 , ^ 2 (* A ) = *i, 



(A.3) 



we obtain the class of 8— self-similar solutions. If we as- 
sume as the simplest case g\ = a, g 2 — A 2 a and 53 = l+5a 
(A and S constants, see Table 1 in VT98 ) we get the solu- 
tions presented in ST94, TT97, STT99 and in the present 
paper. 

Appendix B: Equations for meridional self-similar 
flows) 

From the definitions of Sec. |2.2| we deduce the following 
expressions for the three components of the velocity and 
magnetic field (for details see ST94): 



B r = B, 



1 



G 2 (R) 



cos ( 



Bo = 



B , 



— £>* „„ ; — — - — sin 



-B, 



G 2 (R) 2 

A 1-G 2 (R) 
G 2 (R)1-M 2 (R) 



R sin 9 . 



V r = K 



M 2 (R) 



Vg = -V, 



G 2 (R) y/l + Sa(R,e)' 
M 2 (R)F(R) sint 



G 2 {R) 2 y/i + 6a(R,e) ' 
A G 2 (R)~ M 2 (R) R sin 6 



G 2 {R) 1-M 2 (R) ^l + 6a(R, 



(B.l) 
(B.2) 

(B.3) 
(B.4) 

(B.5) 
(B.6) 



The three ordinary differential equations for H(R), 



M 2 {R) and F(R) are: 



dn 

dR 



2 

"g 1 



dM 2 M 



J2 



dR R? 



(F-2) 



M 2 R 2 ' 



dF(R) _ Af F (R, G, F, M 2 , U; k, S, v. A) 
dR ~ RV(R, G, F, M 2 ; n, A) ' 

Af M (R, G,F,M 2 ,U;k,S, v, A) 



dM 2 (R) 
dR 



RV(R,G 7 F 7 M 2 ;k,\) 



(B.7) 
(B.8) 
(B.9) 
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where we have defined: 

V = (M 2 + R 2 X 2 ^ , (B.10) 



Nf = -{S - k)i 



, RG 2 



2M 2 

(2kHG 2 R 2 + (F+ 1)(F - 2)) x 
R 2 F 2 



F 



M 2 F 

+ ——{F-2)[F + 2 + 2k 



"G 2 



N 2 

2R 2 X 2 ^ 



-X 2 R 2 F(F-2)^§ 



+X 2 R 2 I 1 



I? 

: 'G 2 



N 2 

P 2 \ 2 B 
~Ffi 



D 2 
F 
~2 



N 



2 Ny 



D 



M 2 D 



N m = (5-k)v 2RG2 <** 2 



2M 2 
R 2 



(M - 1) 



^UR 2 G 2 M 2p - - ^-{F - 2)(4k^ + F + 4) 
2 4 G 2 



M 



R 



+ — (F-2)(8«^ + F 



4F 



X 2 R 2 (F - 2) 



D 



' N 2 

X 2 R 2 (2M 2 + F-2) ( 



1 N 2 



D 2 2M 2 D 2 



(B.12) 



with 



1 - G 



M 2 -G 2 , D = l-M 2 . (B.13) 



The meaning of the various parameters is discussed in Sec. 
2. 

At the AlfVen radius, the slope of M 2 (R = 1) is 
p = (2 — F*)/r, where r is a solution of the third degree 
polynomial: 



t 6 +2t z 



kIT* F 2 - 4 



A 2 



4A 2 



- 1 



T+ <^ = 0.(B.14) 



and the star indicates values at R = 1 (for details see 
ST94). 

Appendix C: Forces acting on the outflow on the 
poloidal plane 

The components of the various MHD forces in the poloidal 
plane along the natural unit vectors (n, s) are related to 
those along the spherical coordinates poloidal unit vectors 
(f, 9) via the following expressions: 



f(R, 9) = f cos X (i?, 9) + f sm X (R, 9) , 
f (R, 9) = f sin X (R, 0) ~ f cos X (R, 0) , 



(C.l) 
(C.2) 





Fig. C.l. Sketch of poloidal unit vectors for a converging 
(x < 0) and a flaring (x > 0) streamline (thick line). 



with x the angle a p oloid al streamline makes with the 
radial direction f (Fig. |C.l|) . This angle is related to the 
expansion factor F(R) and the polar angle 9 as, 



tanx(i?, 9) = — tan 



MR) 



(C.3) 



When x < the poloidal streamline turns towards the 
symmetry axis, while when x > it flares faster than ra- 
diality (see Fig. C.l). By taking into account Eqs. (B.l) 
- (B.6) for the dependence of the physical quantities on r 
and 9, we obtain the following expressions for the merid- 
ional components of the various MHD forces f s and /" 
plotted in Figs. || - [l2[ Note that all forces have the di- 
mensions of p*V 2 /2 = B 2 /4ir, a common factor which 
for simplicity we shall omit and the forces will thus be 
given dimensionless in the following (the different forces 
and their direction are listed in Table CI). 
• Gravity: f g = -pQM/R 2 f. 

The two components of the gravitational force are: 



dV 

fg = ~P-q^ = - C0S X(R, a) 
dV 

fg = ~PlT = ~ s[n x(R, a) 
on 



APR 2 



M 2 R 2 



(1 + Sa), 



(1 + 8a) . 



(C.4) 



(C.5) 



Note that always /| < 0, so that gravity in all cases acts to 
decelerate the flow, as expected. Conversely, /" is positive 
and assists further a streamline which is turning toward 
the axis (x < 0), while it is negative assisting further 
decollimation when x > and there is a flaring of the 
streamlines. 

• Pressure gradient force: / VP = — VP(i?,a). 

The two components of the pressure gradient along 
and perpendicular to a poloidal streamline are: 



f 



VP 



dP 



cosx(^, "i^i 1 + Ka ) > 



(C.6) 
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rn 

JVP = 



dP 
dn 



= — sinx(i?, a) 

an 

dR 



4 


dn 


dR ' 


4n 




~FR 


(t 



AkRR 
' FG 2 

- 1 



(C.7) 



Note that f^ P is always positive driving thus the wind 
whenever n decreases monotonically with R, as for ex- 
ample close to the hot base of the outflow (Figs. [)] - |l2|, 
left panels). In the strongly oscillating regions fy p pe- 
riodically accelerates and decelerates the wind. The be- 
haviour of the perpendicular component fy P is more com- 
plex because it depends also on the strength of nil. In the 
transAlfvenic region the pressure gradient is large, so that: 



/v 



VP 



sin x(R> a ) (1 + Ka ) 



dn 

dR 



while for R > 1, dll/dR -> and 
4kR, 



VP 



FG 2 



n sin x(R> a ) 



(C.8) 



(C.9) 



Hence, it becomes evident why in the subAlfvenic (R < 1) 
and transAlfvenic (R >1) regions the gradient of the gas 
pressure acts to 'open' the fieldlines {fy P < 0), while 
further away (R ^> 1) the term cx nil acts to assist colli- 
mation for k > (/y P > 0). 
• Centrifugal 'force': f Vi/> = zb(pV 2 /za). 

This centrifugal 'force' is always directed outwards, in 
the +zb direction with components : 



'2 
'JP 

ZU 



ds 



cosx(R, a) 



2A 2 



RG 2 M 2 



M 2 -G 2 



M 2 



F 

1 

2 



(CIO) 



ril 

JV,„ 



pV£ dw 
zu dn 



sm\{R,a) 



4A 2 



FRG 2 M 2 
G^ + ' 



M 2 - G' 



-i 2 



1-M 2 



F 



- 1 



(C.ll) 



We see that always fy > 0, such that the centrifugal 
'force' acts towards accelerating the flow, as a 'bead slid- 
ing along a rotating wire' in the picture of Blandford & 
Payne (1982). On the other hand, we have fy < such 
that the centrifugal 'force' tends to further decollimate the 
streamlines, for all angles 8 when the poloidal streamlines 
are deflected from radiality towards the polar axis (x < 0, 
or, < F < 2) and also for small polar angles 8 satis- 
fying tan 2 6 < 2/(—F), or, tan 2 \ < {-F)/^ when the 
streamlines flare towards the equator (x > 0, or, F < 0). 

• Toroidal magnetic force, fg = (V x B v ) x B v . 

The toroidal magnetic force fg has two components, 
a tension Tg v and a toroidal magnetic pressure gradient, 
fvB v , both appearing in Eqs. (16) and (|l7|). 

• Magnetic hoop stress of B v : T g v = iSj(B 2 /Attvu). 



The tension of the toroidal magnetic field is in a direc- 
tion opposite to the centrifugal 'force', i.e., it acts in the 
—to direction with components : 



dw 



Attzu ds 



cosx(R,a) 



2A 



2 r 



1-G 2 



RG 2 
1 - 



Tn 
R — 



B v dvu 
4ttvu dn 



— sinx(i?, a) 
R 2 

G 2 ~ + ' 



4A 2 



FRG 2 



1 - M 2 
F 

1-G 2 
1 - M 2 



(C.12) 



(C.13) 



Thus, we see that always T B < 0, such that the hoop 
stress acts towards decelerating the flow. On the other 
hand, we have T p > such that the toroidal magnetic 
stress tends to collimate the streamlines when the poloidal 
streamlines are deflected from radiality towards the polar 
axis (x < 0, or, < F < 2) and also for small polar angles 
9 satisfying tan 2 9 < 2/(—F), or, tan 2 \ < {~F)/2 when 
they flare towards the equator (x > 0, or, F < 0). 
• Magnetic pressure gradient of B v : f^ B = — V(B 2 /8ir). 

This force is positive or negative depending on whether 
Bp decreases or increases with zb, respectively. Its com- 
ponents are: 







2A 2 1 



G 2 



2 \2 



F 

2' 1 



(G 2 



R 



1 - G 2 dM 2 



1-M 2 dR 



M 



a, (C.14) 



Jv 




4A 2 



1-G 2 



FRG 2 (1 - M 2 ) 2 

FR 1 - G 2 dM 2 
~Yl-M 2 dR 



F-G z 



1 



(C.15) 



In the accelerating initial region B v decreases with R and 
thus with zu so the positive gradient of the toroidal mag- 
netic field accelerates the flow. In the case of collimated 
flows B v increases with zu so the negative gradient of the 
toroidal magnetic field acts in the same direction as the 
hoop stress, collimating and decelerating the wind. Note 
that it is easy to see that the meridional components of the 
tension and the gradient of B v have an equal amplitude. 
• The poloidal flow inertial and curvature forces, fy = 
-p(V p -V)V p . 

The two components of this force are: 



- pV *-dS- 



M 2 G 2 
4i? 3 



2 cosx(-R, a) 
G 3 



dM 2 
~d~R 

CP 
R 2 



M^ 
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(F-2) 



F 2 
T 



El 

2 



2F Z - 2F 



dM 2 
~d~R 

. (C.16) 
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Table C.l. Forces acting on a converging (\ < 0) and a flaring (x > 0) streamline: they are parallel or antiparallcl to 
the unit vectors s and n when they are > or < 0, respectively. Then positive or negative f s means that it accelerates 
or decelerates the flow, while positive or negative /" means that it collimates or decollimates. Conversely fy < 
corresponds to acceleration and f v > to deceleration. 
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Components 
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where R c is the local radius of curvature. Along a poloidal 
streamline, f v balances all other forces along s on the 
poloidal plane. Obviously, when the poloidal flow inertial 
force —fy is positive the flow is accelerated, as in the 
inner region, while when it is negative the flow deceler- 
ates, as when oscillations are present. Perpendicularly to 
a poloidal streamline, i.e. along n, the poloidal flow cur- 
vature force —fy is positive when the streamlines turn 
towards the axis (as in the inner region), zero in the ra- 
dial case and negative when there is some flaring away 
from radiality (e.g., in the oscillating zone). 

• Poloidal magnetic force, f b p = (V x B p ) x B p . 

This force is directed only normal to the poloidal mag- 
netic ficldlincs, since its s— component does not appear in 
Eq. (16). As for B v , the poloidal magnetic force has two 
components, a tension Tg^ and a poloidal magnetic pres- 
sure gradient, /vb p > both appearing in Eq. (p"7|). 

• Magnetic hoop stress of B p : Tg = n ■ (B p ■ \?)B p /4ir. 

This force is focusing the streamlines towards the axis 
when they collimate and towards the equatorial plane 
when they are flaring: 



rpn 
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El 
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'F*_ dF_" 
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(C.18) 



Note that in the regime of the oscillations this poloidal 
magnetic curvature force is in opposite phase with the 
poloidal velocity curvature force. However, the ratio of 
their amplitude equals to M 2 and thus in highly super- 
Alfvenic outflows (M ^> 1), the poloidal magnetic curva- 
ture force is negligible in comparison to the poloidal flow 
curvature force. 

• Magnetic pressure grad. of B p : f^ B = — ri-V (-B p 2 /87r). 

The pressure gradient of the poloidal magnetic field 
reads as: 



2(2 - F) 




(C.19) 



In an almost radial expansion B p drops rather fast with 
R, B p r~j R~ 2 and thus f^ B has a strong component 
along the radial direction. In such cases, the gradient of 
the poloidal magnetic pressure has a component along - 
n acting to decollimate the outflow. In fact, this strong 
poloidal magnetic pressure gradient when combined with 
the rather weak poloidal magnetic curvature force for su- 
per Alfvcnic flows, results in a total poloidal magnetic force 
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which has roughly the magnitude and sign of /y B , Figs. 

re 



